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Abstract 

For square contingency tables with ordered categories, the present paper gives the theorem that the symmetry model holds if 
and only if four models of (1) cumulative sub-asymmetry, (2) palindromic symmetry, (3) global symmetry, and (4) marginal 
means equality, hold simultaneously. It analyzes some data using this decomposition. 
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Introduction 


Consider an rxr square contingency table with ordered categories. Let P u denote the probability that an 
observation will fall in the i th row and j th column of the table (i = !,■■■, r;j = !,■■■, r). Consider the symmetry (S) 
model, defined by 

Pa = Pa O' * i\ 

(Bowker, 1948; Bishop, Fienberg and Holland, 1975, p.282). 

Let 

G ,j = Z Z P* 0 < j), 

5= 1 t- j 

and 

G Ji = Z Z Pa 0 < j). 

5=1 t- j 


The S model can be expressed as 


G,=G, (/<./>. 


McCullagh (1978) considered the palindromic symmetry (PS) model, defined by 


G, - A tX ‘ G, (1 <i<j<r), 
a i-i 

where a, = 1 without loss of generality. A special case of the PS model obtained by setting A = 1 and «, = • • ■ = a r , 
is the S model. 

Let X and Y denote the row and the column variables, respectively. Consider the marginal means equality (ME) 
model, defined by E(X ) = E(Y) , where 


E(X) - Z Z 

i=l j = i 


'P 
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and 


Also, the ME model can be expressed as 


E(Y) = ±±jP, r 

<=i j = i 


Lg„ +1 = £g, h 


(see, Tahata, Yamamoto and Tomizawa, 2012). 

As a model which indicates the structure of {g, 7 } with |' j -i\= 2 , namely, j/r.j with j - i |> 2 , Tomizawa, 
Miyamoto and Ouchi (2006) proposed the cumulative sub-symmetry (CSS) model, defined by 

G iji+2 = G i+%i (i = 1, 2, - - ■ , r — 2). 


Tahata et al. (2012) gave the theorem that the S model holds if and only if all the PS, ME and CSS models hold. 

Tomizawa, Miyamoto, Yamamoto and Sugiyama (2007) considered the cumulative two ratios-parameter symmetry 
(C2RPS) model, defined by 

-r©' ■ </<./>. 

G n 

The global symmetry (GS) model is defined by 

= s L , 

where 

S u=YY.P<p S l = YLPp 

i<j i<i 

(Read, 1977). Tahata, Yamamoto and Tomizawa (2013) gave the theorem that the S model holds if and only if all 
the C2RPS, GS and ME models hold. 

The marginal homogeneity (MH) model is defined by 

G„ +1 =G, +U 0 = 1, 2, • • • , r - 1 ), 


(Stuart, 1955). The model indicates that the marginal distribution for the row variable is equal to that for the 
column variable. Note that the model indicates the structure of {G^.j with | j-i |= 1 . The extended MH model is 
defined by 

G„ +l = 8G M l (i ~ 1 , 2, - - - , r — 1 ), 

(Tomizawa, 1993). This model indicates the asymmetric structure of j G tj | with \ j -i\= \ . Therefore, we are 
interested in considering the extended CSS model as the analogous to Tomizawa (1993). In this paper, we propose 
a model which indicates the structure of asymmetry for j G /( j with | j-i |= 2 . Also, we show that the C2RPS model 

can be separated into the proposed model and the PS model. Using this result, we give the decomposition of the S 
model. 


Decompositions 

This section proposes a model based on jG^.j , and gives a decomposition of the C2RPS model. 


Consider a model defined by 


= H (« = 1, 2, • • ■ , r - 2). 
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We shall refer the model as the cumulative sub-asymmetry (CSAS) model. We now obtain the following theorem. 
Theorem 1. The C2RPS model holds if and only if both the PS and CSAS models hold. 

Proof. If the C2RPS model holds, then the PS model holds. Assume that the C2RPS model holds, and then we shall 
show that the CSAS model holds. We see from the C2RPS model that 

= r©" +2 >-' = r © 2 0 = 1, 2, — ,r -2). 

G i+2,i 


Thus the CSAS model holds. 

Assume that both the PS and CSAS models hold, and then we shall show that the C2RPS model holds. We see from 
the PS model that 


Wi +2 A a i r 1 O 
— = A 0 = 1, 2, — 2). 

G i+2,i a i+\ 


Then, we can set a, = 1 without loss of generality. Since the CSAS model holds, 

— h- = a — = A — = Q. 


So, we obtain 


Similarly, we obtain a 3 = (A/Q) 2 because 


Therefore, we can see 


a, = — . 

- Q 


^21 = A — = = Q. 

G 42 cl 3 Q a 3 


01 M ~ Q (* = 1, 2, — , r — 2). 


Thus 



n 



a < j). 


The C2RPS model holds. The proof is completed. 

We now obtain the following theorem from the result given by Tahata et al. (2013) and Theorem 1. 
Theorem 2. The S model holds if and only if all the PS, CSAS, GS and ME models hold. 


Goodness-of-Fit Test 

Let h,. denote the observed frequency in the i th row and j th column of the rxr table with n - /, T" n n , and let 
ntij denote the corresponding expected frequency. Assume that a multinomial distribution applies to the rxr table. 
We can test each model for goodness-of-fit by using the likelihood ratio chi-squared statistic G 2 defined by 


° 2 = 2 Z Z ,! y lo 8 

i=l j=l 




where my is the maximum likelihood estimate of under the model. Table 1 presents the numbers of degrees of 
freedom for testing goodness-of-fit of the models. 
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TABLE 1 The DEGREES OF FREEDOM OF MODELS FOR r x r TABLE 


Models 

Degrees of freedom 

S 

r(r — 1)/2 

C2RPS 

<N 

1 

<N 

7 

PS 

CN 

1 

Sn. 

7 

CSAS 

r-3 

ME 

1 

GS 

1 


An Example 

Table 2, taken directly from Tominaga (1979, p.130), is the cross-classification of Japanese father's and his son's 
academic background which were examined in 1955. Note that category (1) is elementary school; (2) junior high 
school; (3) high school; and (4) university. Table 3 gives the values of the likelihood ratio chi-square statistic G : for 
models applied to these data. The CSAS model fits the data in Table 2 well. 


table 2 Cross-classification of Japanese father's and 

HIS SON'S ACADEMIC BACKGROUND; TAKEN FROM TOMINAGA (1979, P.130) 


Father's 

Son’s educational level 

Total 

educational level 

(1) 

(2) 

(3) 

(4) 

(1) 

374 

a (374.0000) 

602 

(602.0000) 

170 

(167.3926) 

64 

(64.6918) 

1210 

(2) 

18 

(18.0000) 

255 

(255.0000) 

139 

(139.0000) 

71 

(72.9156) 

483 

(3) 

4 

(6.9571) 

23 

(23.0000) 

42 

(42.0000) 

55 

(55.0000) 

124 

(4) 

2 

(1.5482) 

6 

(3.4947) 

17 

(17.0000) 

53 

(53.0000) 

78 

Total 

398 

886 

368 

243 

1895 


Note: a Estimated expected frequencies for the CSAS model. 


Under the CSAS model, the maximum likelihood estimate of Q is Q = 27.2872 . Hence, under this model, the 
probability that a father's academic background status is i or below and his son's academic background status is 
i + 2 or above is estimated to be Q = 27.2872 times higher than the probability that the father's academic 
background status is i + 2 or above and his son's academic background status is i or below ( i = 1,2) . 

Moreover, from Table 3, we see that the S model fits these data poorly, yielding G 2 =1151.23 with 6 degrees of 
freedom. Also, each of the PS, ME and GS models fit these data poorly, but the CSAS model fits these data well. By 
using the decomposition of the S model (i.e.. Theorem 2), we shall consider the reason why the S model fits these 
data poorly. We can see that the poor fit of the S model is caused by the poor fit of the PS, ME and GS models 
rather than the CSAS model. 


table 3 Likelihood ratio chi-square values G 2 for models applied to the data in Table 2 


Models 

Degrees of freedom 

G 2 

S 

6 

1151.23* 

C2RPS 

4 

31.20* 

PS 

3 

31.13* 

CSAS 

1 

3.18 

ME 

1 

700.17* 

GS 

1 

1093.22* 


Note: *Significant at the 0.05 level. 
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Concluding Remarks 

We have proposed the CSAS model which indicates the structure of | G, y ] with j -i\=2 , namely, j p tj j with 
| j - i |> 2 . Using this model, we have given Theorem 1. The theorem may be useful to see the reason for the poor fit 
when the C2RPS model fits the data poorly. In practice, we can see from Theorem 1 that the poor fit of the C2RPS 
model is caused by the poor fit of the PS model rather than the CSAS model in example (see Table 3). For square 
contingency tables with ordered categories. Theorem 2 would be useful for seeing the reason for the poor fit when 
the S model fits the data poorly. 
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